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Abstract 

Unitary 1-matrix models are shown to be exactly equivalent to hermitian 1-matrix models 
coupled to 2N vectors with appropriate potentials, to all orders in the expansion. This 
fact allows us to use all the techniques developed and results obtained in hermitian 1-matrix 
models to investigate unitary as well as other 1-matrix models with the Haar measure on 
the unitary group. We demonstrate the use of this duality in various examples, including: 
(1) an explicit confirmation that the unitary matrix formulation of the M = 2 pure SU{2) 
gauge theory correctly reproduces the genus- 1 topological string amplitude (2) derivations 
of the special geometry relations in unitary as well as the Chern-Simons matrix models. 
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1 Introduction 



Random matrix models P have served as cornerstones for various areas of theoretical physics. 
Hermitian matrix models, particularly those in double-scaling limits, were intensively studied 
in the early 90's in connection to two-dimensional quantum gravity |2j. Unitary models, on 
the other hand, are known to arise as a one-plaquette lattice model of two-dimensional 
Yang- Mills theory More recently, matrix models, both hermitian and unitary ones, have 
proven to be powerful tools for studying supersymmetric gauge theory and topological string 
theory [HE]. 

Both hermitian and unitary 1-matrix models are reduced, after the gauge fixing, to eigenvalue 
models with different Jacobians. Despite the difference, they induce at short distances the 
same interaction between the eigenvalues, and share many common properties. For instance, 
the critical behaviors of both unitary and (even-potential) hermitian 1-matrix models are 
described by the Painleve II equation when two clusters of eigenvalues get coalesced jB] , and 
more generally they are governed by the mKdV hierarchy [7]. 

The recent use of matrix models for the study of gauge theory and string theory [HE] requires 
not only the knowledge of their critical behaviours but also their individual higher genus 
corrections away from criticality; the technology to compute them has been less developed 
in unitary 1-matrix models than in hermitian ones. Moreover, some unconventional matrix 
models, which also have the Haar measure on the unitary group, have recently been proposed 
to describe a Chern-Simons theory [HI El and five-dimensional gauge theories [TUI. Therefore 
it is useful to develop computational techniques in these models as well. 

In this paper we show that, in fact, unitary and hermitian 1-matrix models are nothing but 
different descriptions of exactly the same system^ ; by a simple change of variables we map a 
general unitary 1-matrix model to a certain hermitian 1-matrix model coupled to 2N vectors 
with an appropriate potential. A similar transformation can also be done on the Chern- 
Simons matrix model. This means that all the techniques developed and results obtained 
in hermitian 1-matrix models — the higher genus computations of partition functions and 
correlation functions [121 CSl lEl US] , the relations to integrable hierarchies [THIITTIITHIIT^I^ ^ 
and conformal field theories pi] , the derivation of the special geometry relation (221 — 
they all can be used to investigate unitary 1-matrix models as well as other analogous models 

^In the double scaling limit, the duality of the two matrix models was established more than a decade 
ago. The main point of the present paper is to show that the duality holds even off criticality so that one 
can use it to compute higher genus topological string amplitudes, etc. 

^For unitary 1-matrix models it has been known that their underlying integrable structure is the "modified 
Volterra" [TQ or the "quarternionic Toda" ^Jj hierarchy. 



1 



with the unitary measure. 



The relevant transformation is simply the fractional linear transformation 

- - <-) 

where U and $ are a unitary and a hermitian matrix, respectively. The fact that their Haar 
measures are related by 

may be found in the book [24J originally published in 1958. This transformation was used 
in to derive the Virasoro constraints |2S] in unitary 1-matrix models, although the fact 
that it implied an exact duality does not seem to have been fully recognized ^ . In pH] the 
possibility of using the transformation (jl.lj) was mentioned in the study of matrix models 
coupled to an external field, with a conclusion that it did not have "a particularly attractive 
expansion" for their purposes ^. Later on, the techniques of computing the higher genus 
corrections were developed in hermitian 1-matrix models fT^ 1131 ITU I15j . and more recently 
were used in the Dijkgraaf-Vafa theory However, the transformation ()1.1|) has never 

been utilized in this context. 

The dual hermitian matrix models we end up with are closely related to what we encountered 
in the study of four-dimensional Af = 1 U{n) gauge theories coupled to Nf fundamental 
matter fields (e.g. 1231221; see also j2S| for further references); in the present case we will 
have vectors of order (= 2A^), however. Nevertheless, we can take advantage of the 
derivation of the special geometry relation given in [23] to prove analogous relations for 
unitary and the Chern-Simons matrix models. 

The organization of this paper is as follows: In Section 2 we describe the details of the 
unitary/hermitian duality transformation, which plays a central role in the present paper, 
and demonstrate the use if it in the 1-cut Gross- Witten model as an illustration. In Section 
3 we establish the relationship between the resolvent commonly used in unitary 1-matrix 
models and that in the dual hermitian models. We then derive the special geometry relation 
for general unitary 1-matrix models, following j2Sl- Section 4 and 5 are devoted to the 
applications of the duality transformation. In Section 4 we use it in the unitary matrix 
formulation of the M = 2 pure SU (2) gauge theory and prove that it correctly reproduces 

■^It was noted in the same paper f^l9j. Appendix D) that the question of whether the partition function 
of the full unitary model was a r function of any integrable hierarchy remained open. 

'''In there was also mentioned another possibility U = e'* [2ZI , which was recently exploited by Tierz 
P5| in the study of the Chern-Simons matrix model. See section 5 for more details. 
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the genus-1 topological string amplitude, whereas in Section 5 we deal with the matrix 
model mirror of Chern-Simons theory proposed in I8j. Finally we conclude the paper with a 
summary. The two Appendices contain some technical details needed in the text. 



2 Unitary/Hermitian Duality 



2.1 The Duality Transformation 



The partition function of the unitary 1-matrix model with a potential Wu{U) is 



1 



/ 



(2.1) 



Z 



Vo\U{N) 



U is an N X N unitary matrix and dU is the Haar measure on U{N). is the coupling 
constant. It is a classic result that this matrix integral is reduced to that over the eigenvalues 
120] as 



where we take — vr < < vr (j = 1, . . . , N). 

The partition function of the hermitian 1-matrix model with a potential iy($) is given by 



$ is an X hermitian matrix and Xj (j = 1, . . . , A^) are the eigenvalues of 

In hermitian models the Jacobian associated with the gauge fixing is known to be a Vander- 
monde determinant; the log of it may be viewed as a potential for a repulsive force between 
two eigenvalues on the real axis. For unitary models, on the other hand, sin replaces the 
linear function in the Jacobian, so that it may also be viewed as representing the same linear 
force between two eigenvalues of U, located this time on a unit circle because 



Z 




(2.2) 




(2.3) 



(2.4) 



. 2^k - 

m 

2 



(2.5) 
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To write Z as ^hermitian for some we have only to use the fractional linear transforma- 

tion 

i — z , ^ 

w = 2.6 

t + z ^ ' 

on the complex plane, which maps the unit circle \w\ = 1 to the real axis z G R. U w = e*^, 
it is equivalent to 

z = tan-. (2.7) 
With this elementary change of variables the determinant factor becomes 

N a _ o N 1 N 

n-n^\^ = n7T-4^n(A.-Ao^. (2.^ 

k<i ^ j=iU + ^jJ k<i 

Since 



ie, = ^ (2.9) 



has another factor of ^t)?, we find 



z = 2"/n«,e-="-"(Si)-"'-"**''n(A.-A, 



2 

I) ■ 



(2.10) 



j=l k<l 

Therefore, in the 't Hooft expansion [3T| with 

/i = g,N (2.11) 

kept fixed, factors out in the exponential, leaving a hermitian matrix model with the 
potential 

W{X) = Wu('-^)+i^log{l + X'). (2.12) 



Z ()2.10|) may be thought of as arising from a hermitian matrix model coupled to vectors: 

N 



1 f ' 

^+vector - Yo\U{N) 7 K ^Qi^QI exp 



1=1 



1=1 
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oc |rf$e-^(""^*^"^"^^^-"^^"-'"^0 (2.13) 
after integrating out Qj, Q\. Since 

/xlog(l + A2) = ^,ivflog(A + + log(A-0V (2.14) 



we see that the unitary 1-matrix model partition function Z ()2.2|) is, up to an A^-independent 
multiphcative constant, equal to that of the hermitian matrix model Z+vector ()2.13p with a 
potential 



coupled to N I = 2N vectors having pure imaginary "masses" ±i. 
Some remarks are in order here: 

1. The potential W{\) (j2.12j) does not depend on N. Therefore the equivalence exactly 
holds to all orders in the expansion. 

2. The extra log term in W{X) can be expanded as a Taylor series near A = 0. Thus all 
the techniques developed for the computations of higher genus amplitudes O HSI HI] 
can also be applied ^ to those of unitary 1-matrix models. This is in contrast to the 
transformation of the type such as U = expi$, which leads to a logarithmic singularity 
at A = 0. 

3. An immediate consequence of the duality is that the partition function of a unitary 
1-matrix model is also a r-function of the KP hierarchy. 

4. Z+vector ()2.13p is the very matrix model used for computing the effective superpotential 
and other F-terms of a four- dimensional JV = 1 U{n) gauge theory coupled to Nj 
fundamental matter fields (e.g. [201 122]), although A'^^ = 2A^ > 2n is the region where 
this interpretation breaks down. 

2.2 Example: The Gross-Witten Model 

To demonstrate the use of the duality transformation we will first rederive the well-known 
1-cut solution of the Gross-Witten model j3] by using standard hermitian matrix model 
techniques j^. The Gross-Witten model is characterized by the potential 




(2.15) 





(2.16) 



The corresponding potential of the hermitian matrix model is 



1 + 2;^ 



(2.17) 



^as far as Wu{U) is regular at J7 = 1. For example, Wu{U) can be any polynomial of U and U 
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The resolvent is given by 



_ jiE^}Eli,, ;^^'<-) (2.18) 

47r^/i J A _ a;) y (a; _ «) - 6) 
where a, 6 are the end points of the cut satisfying 

, = 0, 

J {x — a){x — h) 

dx , = 47^^/i. (2.19) 

^ Y'(x-a)(x-6) 

The contour A surrounds the cut, running close enough to it so that the contour avoids the 
branch cuts of log(l + z^). The conditions ()2.19|) follow from the large \z\ behavior of uj{z). 
They are easily solved to give 



^ '> 0), (2.20) 



2-/i 



so that 



u{z) = ^{w\z)-^3^^^\. (2.21) 

The spectral density p(A) may be read off from the discontinuity of the resolvent uj{z) as 

p(\)dX = — —(uj(X + iO)-uj(X-iO))dX 
2Txi 



4 y/fe^ _ ^2 

vr/i v/rTF(l + A2)2 



= Acos-W^-sin^-rf^, (2.22) 
vr/i 2 V 2 2 ' ^ ^ 

which is precisely the solution of the Gross- Witten model in the weak-coupling regime [S]. 
The genus-0 free energy Fq can be obtained as that of the hermitian matrix model |i32j : 

Fo = r dXp{X)W{X)+ I? r rfAp(A)log|A| - ^W{<d). (2.23) 

2 J — OO J —CO 2 

The genus-1 free energy Fi is also known to be given by [12| 

Fi = -^(logMi + logJi + 41og26) (2.24) 
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with 

dw ^W'{w) 



Ml 



A 2711 ^vfi — h'^{w + h) ' 

r dw -W'iw) 
Jl = i — , ^ . (2.25) 



Ta 2m y/w"^ - h^iw - h) 
Straightforward residue computations yield 



^ fA2m -b^(w + b) ^ ■ ' 

4 

Jl. (2.27) 



Therefore 



Fi = ^log(l + 62)-llog2 

= llog(2-/i)-^log2, (2.28) 



which also agrees with the known result of [211 ^ constant. 

3 Resolvents and The Special Geometry Relation 

3.1 Resolvents in Unitary/Hermitian Matrix Models 

In this subsection we will examine the difference between the two resolvents: one is that 
commonly used in unitary 1-matrix models, and the other is the usual resolvent in hermitian 
models. The saddle point equation derived from (j2.2|) is 

1 ^ a — a 

l^cot^ = fe^^W^e^^O- (3.1) 
The resolvent uju{x) is customarily defined in unitary matrix models as 

uu{x) ^ ^Ecot^ ''=T r depu{e)cot^. (3.2) 



Nf^^ 2 ' 2 

On the other hand, the resolvent in hermitian 1-matrix models is given by 
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The relation between the two is as follows: The spectral densities satisfy 

dOpuie) = dXp{X). (3.4) 

It is easy to show that as a differential 

x-e , / 2 2z \ , , 

dxcot^ = dz(^j—^-^^), (3.5) 

where z = tan | and A = tan |. Thus we obtain the relation 

LJu{x)dx = (2uj {z) - {plog{l + z^)y)dz. (3.6) 



If we define 



then we find 



yuix) = -^Wuie''') - puuix), 
dx 

y{z) = W\z)-2pu{z), (3.7) 



yu{x)dx = y{z)dz. (3.8) 



Therefore, the function yu{x) (or rather the differential yu{x)dx), which characterizes the 
curve of the unitary matrix model, is simply computed as y{z) (or y{z)dz) in the correspond- 
ing hermitian matrix model. 



3.2 The Special Geometry Relation 

One of the properties of matrix models that makes contact with gauge theory and topological 
string theory is the "special geometry" relation a relation between the genus-0 free energy 
Fq and periods of the differential y{z)dz. It is due to this property that makes possible to 
identify Fq as a prepotential of special geometry [22] • 

This relation was recognized in For hermitian 1-matrix models with a polynomial po- 
tential, a proof using a Legendre transformation was given in [22] and a more direct one 
in [22] • Also for unitary 1-matrix models, one may argue that a similar relation holds [3] 
because y{x) is again interpreted as a "force" on an eigenvalue exerted by other ones located 
on the unit circle. Strictly speaking however, the relevant period integrals diverge without a 
regularization, the detail of which is not determined by the intuitive argument. In this sub- 
section we will derive, following (and appropriately modifying) the proof in , the special 
geometry relation for unitary 1-matrix models. 
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Let Wu{U) be any polynomial of U and U ^. The genus-0 free energy for n-cut solutions is 
given by 

11 POO " roo 1 

L J -co J~oo I 

with W{X) given by (ITT^ . and 

fij = dXp{X) (j = l,...,n). (3.10) 

a2j-i, a2j are the two end points of the jth cut, whereas Xoj is an arbitrary point on the jth 
cut, for j = 1, . . . ,n. We assume that they are all real and ai < a2 < ■ ■ ■ < a2n-i < CL2n- This 
expression for the free energy is more or less standard, but we give a derivation in Appendix 
A for completeness. 

The right hand side of (j3.9p can be written in terms of contour integrals as 

^0 = ^E/ dzy{z)W{z) + \im f.^ i dzy{z) - ^iy(Ao) + /.^ log |Ao 

(3.11) 

The definitions of the contours are shown in Figure 1. 

y{z) is given by the equation 

y[z) = W'{z) - 2puj{z) 



SO that 

dzy{z) = - Anipj (j = 1, . . . , n). (3.13) 



Our aim is to express |^ in terms of the periods of y{z)dz. For this purpose we need to 
compute ^M^dz, which may be found by simply studying the locations of the singularities, as 
was done in j2Sl- In the present case, however, W'{z) has extra pole singularities ai z = ±i 
with residue p coming from the unitary Haar measure and must be also taken into account. 

y{z)dz has the following properties: 



• Since the contours of ()3.12|) do not enclose x = z, the integrals give single-valued 
functions of z on the complex plane. Therefore y{z) reverses its sign under the exchange 
of the first and second sheets. 
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Figure 1: The contours. Those on the second sheet are shown in broken hues. 



• Since we have assumed that Wu{U) is a polynomial of U and U~^, y{z)dz is regular 
at 2r = cxD, 6b, where we denote here and henceforth a point on the second sheet by its 
value of z with tilde. 

• Since uj{z)dz must be regular at z = ±i on the first sheet, y{z)dz exhibits the same 
singular behavior as W'{z)dz at z = ±z, and as —W'{z)dz at z = ±z. 

These properties, together with the condition ()3.13|) . uniquely determine the /Xj derivative 
of y{z)dz as 



dyjz) 



dz 



Airi LJj{z), 



(3.14) 



where 



(j = l,...,n-l) 



'n 




(3.15) 



tUj's are holomorphic differentials normalized as 




(3.16) 
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and uja--b{z) denotes a zero A-period abelian differential of the third kind having simple poles 
at z = a,b with residue ±1, respectively. 

Having computed ^^j^dz, we are now able to derive the special geometry relation. Differen- 
tiating the free energy Fq by fij, we have. 



diii 



-^Ef ^^WW^W + ^E/ dzyiz)\ogiz' + l) 

" /■ If 1 

-mJ2f^j L uJi{z) + - I dzy{z) - -l^(Ao) 



+ lim 

Aq— »— oo 



/i 



log(A^ + l) +2/ilog|Ac 



-9 E f uJi{z)W{z) + lim -vrz^/ij / 



uAz] 



+ lim 



(3.17) 



This expression can be simplified because one can show that the following relation holds: 



n „ n „ 



'^AZ\ 



lim 

Ao— »— 00 



1 



2 — e p~i-~e 

lim ( / + / + 



- - e) - W{-i - e) 



(3.18) 



4 e^O VJb, Ao Ao 

This relation may be derived in the way parallel to that in by using Riemann's bilinear 
relation; we give a brief sketch of it in Appendix B. Thus we obtain the final result 

5Fo 1 



lim 

2 e^o 

- lim 

2 e-*0 



I — e f — i—e 
+ I + 



lim 



dzy{z) -W{i- e) -W{-t- e) 



jith cut) dzy{z) — W{i — ie) — W{—i + ie) 



(3.19) 



where in the last line the integration contour passes through the ith cut. If we go back to the 
ty-plane ()2.6|) . the contour runs from a point infinitely close to = to another io w = 00. 



4 Application I : = 2 Gauge Theory via Unitary 
Matrix Models 

It is now well-known that the effective superpotential and other F-terms in A/" = 1 super- 
symmetric gauge theory can be calculated using hermitian matrix models IH E] . The M = 2 
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limit may also be achieved by turning off the superpotential with the effective superpoten- 
tial kept extremized with respect to the glueball fields jSl]. On the other hand, it has been 
shown [S] that there is another route to realize more directly the pure SU (2) Seiberg-Witten 
curve via a unitary matrix model without referring to any A/" = 1 theory. In this section 
we will apply the duality transformation to this unitary matrix model to examine whether 
the genus-1 topological amplitude, obtained by other means, is correctly reproduced in this 
framework. 



4.1 2-cut Solutions of The Gross- Witten Model 

To realize the pure SU (2) Seiberg-Witten curve we again consider the Gross- Witten model, 
though this time focus on 2-cut solutions. In this case it is convenient to shift 6 by 90 degrees 
in ()2.16|) and consider the potential 



W{z) = - + /ilog(l + z2). (4.2) 

I I 



Wu{e'^) = -2sin0, (4.1) 

Az 

We let ai < ■ ■ ■ < a4 be the end points of the cuts; they are constrained by 

+ i ] dx I ^ ' = 0, 



A2 



W'{x) 






xW'{x) 




\/nLi(^ - 


ak) 


x^W'ix) 




M=i{x - 


ak) 



dx , = 0, (4.3) 



+ f dx . I = Aniu. 

Ai JAzJ ' - ' 

One can straightforwardly compute the differential y{z)dz to find 



2TII \J Aaa Jx=i Jx = ~-iJ I 



z - x)JY{k=i{x - ak) 



8 \/nLi(^-Qfc) 

C {Z^ + 1)2 ' 



where we put \l]\k=i{i - ak) = JUk=i{~'^ - ak) = c. 



(4.4) 



To incorporate the end-point conditions ()4.3p we go back to the ly-plane ()2.6|) . The two cuts 
are mapped to two arcs located on the unit circle with end points given by 

T. = ^ (A:=l,...,4). (4.5) 
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Then the conditions ()4.3p are fulfilled if and only if T^'s are the roots of the following equation 
of U: 



where 



- 2ifiU'^ + a2U^ + 2iiiU + 1 = 0, 



(72 = T1T2 + T1T3 + T1T4 + T2T3 + T2T4 + T3T4 



(4.6) 



(4.7) 



is a free parameter, which is eventually identified as (essentially) the u modulus of the 
Seiberg-Witten curve. Using the relations between the roots and the coefficients, y{z) ()4.4|) 
may be written in the form obtained in 0: 



y{z)dz = 2 J cos^ 9 + /i(sin 9 — u) dO 



with u = The AT = 2 curve and differential are obtained by restricting the moduli to 

a co-dimension one sublocus in the following way: We first scale y{z)dz as 

y{z)dz = 2^ysw{0)de (4.9) 

and take the /i — > 00 limit. Then we obtain the Seiberg-Witten curve for the pure SU{2) 
gauge theory: 



yl^ -sinO + u = 0. 



(4.10) 



On the other hand, we choose, as was done in [HI], Hi and fi{= + (12) (instead of fii and 
fi2) as independent parameters. Using the special geometry relation in the previous section, 
we have 



dFo 




dFo 




dFo 










M2 


dii2 


Ml 



where B = 52-8^^. Therefore yswi(^)dO is indeed a Seiberg-Witten differential. 



(4.11) 



4.2 Gravitational F-terms via The Unitary Matrix Model 



We are now interested in the genus-1 free energy of this matrix model. In the 2-cut case it 
is given by 



F^ = log f[Mj-^ log \K{k)\ - ^ log n(ai - a^)' + \ logd^i - ^sl 1^2 - 04I), (4.12) 



13 



where the moments are defined by 

2m \Ja, JaJ {x-a,)Jut=i{x-ak 



Mi = T^.[j> + i ) dx ^ I (4.13) 



and 



t2)(l - A;2t2) 

is the complete elhptic integral of the first kind with modulus 



Kik) = dt , ^ (4.14) 

^0 Jil-mH-kH^) 



^2 ^ (qi - 04) (a2 - Qs) ^ 2Vm2 _ X 

{ai - an){a2 - tti) u + \/v? - 1 ' 

Here we note that ?/|^ ()4.10|1 has at most only two zeros on the unit circle, and hence is only 
defined as an analytic continuation of the 2-cut solutions. Therefore we evaluate Mj's where 
it is well-defined, and then take the /i ^ 00 limit afterwards. This is very different from 
the situation in the hermitian matrix realization [33] where a spectral curve "really" exists 
in the J\f = 2 limit. Nevertheless, despite the difference, Fi ()4.12|) obtained as an analytic 
continuation correctly reproduces the known result, as we will see below. 

Ml is found after some straightforward residue computations as 

2 r 

Ml = r— "2 — — - — 8 + 6(02(13 + 0204 + 0304) — 3ai(a2 + 03 + 04 — 02^304) 



Hc^{ai + 1) 



2(-/iai + 2) 



+ai(-5 + 3(0203 + a20-4 + a3Ct4)) — 77- (4.16) 

J fj,c{ai + 1) 

with = Various symmetric polynomials of a^'s turn out to be 

2/i 

ai + a2 + 03 + 04 = -, 

fju — 1 

4 

aia2 + OiOs + aitti + 0203 + 0204 + a^a^ = 2 + 



010203 + 01O2O4 + 01O3O4 + O2O3O4 



/iU — 1 ' 
2/i 



jju — 1 

01O2O3O4 = 1. (4-17) 



Thus we find 



n = = /iV (1 + Oif^-')) . (4.18) 

j=i 
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The discriminant can be computed as 

16V'(/i'-4/iu + 4)2(^2-1) 



Y{{ai - a 



i<j 

162(1x2 - 1 



{fXU - 1)6 



/i— »oo 



(4.19) 



and also 



|ai - agl |a2 - 04! = 2 ( 1 + y 1 - — j . (4.20) 

Plugging these data into ()4.12|) . we see that the logw terms precisely cancel. This leaves 

Fi = -^logK(A;)-^log(M2-l) + ilog(M + v/^I^)-^log2/x + 0(/x-^), 

(4.21) 

the M-dependent part of which is in agreement with the known genus-1 result ISEIISSI- 



5 Application II : A Matrix Model Mirror of A Chern- 
Simons Theory 

5.1 The Transformation 

The final example we consider is a matrix model mirror of a Chern-Simons theory. It has 
been proposed that the partition function of the SU{N) Chern-Simons theory of is 
equal to the partition function of a matrix model with a "Wick-rotated" Jacobian as 

Zcs = e-^(^^-^)Z, 

Z = l[du,e-t^Y[smh^^^^^-^. (5.1) 

•' j=i k<i 2 

In this section we apply the duality transformation to this model, compute the genus-0 and - 
1 free energies and confirm that they indeed reproduce those of the Chern-Simons theory. 
In fact, it was already proven by Tierz |2S] that the partition function ()5.ip is equal to the 
Chern-Simons partition function to all orders in the 1/N expansion, by utilizing yet another 
transformation (as in footnote to map ()5.H) to that of a hermitian matrix model (different 
from ours) and use its orthogonal polynomials. However, there are still some benefits in 
considering our transformation in this model as follows: First, it will also demonstrate how 
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to compute explicitly the lower genus free energies in other such models that have the same 
measure of the "Wick-rotated" type, such as those proposed to be related to some five- 
dimensional gauge theories |lOj, even where no good orthogonal polynomials are available. 
Secondly, because our transformation induces no singularity near the origin in the potential, 
one can show the special geometry relation in these models by only slightly modifying the 
proof for hermitian models. Finally, for the same reason, the relation to the KP hierarchy is 
made manifest by our transformation, and not by the one adopted in j2Hl because of its log 
branch at the origin. 

For the measure of the type ()5.ip the transformation we need is 

Xj = tanh^. (5.2) 

Z (15.111 can then be written as 



N 



N 



Jn(A.-AO^. (5.3) 

j=l k<l 

Comparing it with ()2.13j) . we find that the Chern-Simons matrix model has, up to an A^- 
independent multiplicative constant, the same partition function as that of the hermitian 
matrix model with a potential 

W{z) = \ (log[^)' + /xlog(l - z^), (5.4) 

or as that of the hermitian matrix model coupled to N f = 2N vectors having real masses 
±1 with a potential 

2 V"" 1 - 



W^iz) = l-(\ogl^) . (5.5) 



In this case, therefore, the potential of the resulting dual hermitian matrix model has branch 
cuts of log ending on the real axis; this is not a problem because the eigenvalues are located 
only on the interval between 1 and —1, as is clear from the transformation ()5.2|) . 



5.2 The Spectral Curve 

We will now rederive the (nontrivial part of the) geometry of the deformed conifold as the 
spectral curve, following the standard procedure as before. Let us consider a 1-cut solution 
with a cut [a,b]. The end point conditions ()2.19p with W ()5.4|) imply 

b = -a = Vl- e-/^ (>0). (5.6) 
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The resolvent is also given by ()2.18p with W ()5.4|) . The contribution of the second term of 
W is simply obtained by residue computations, while the contribution of the first term can 
be written as 



where 



Kz) 



dx- 



21og(l + a;) 



(5.7) 



(5.8) 



If the contour A is deformed to a large circle around infinity, I{z) also picks up the contri- 
bution from the difference of the integrals along both sides of the branch cut of log(l + x). 
After some calculations we find 

Am ^ z + b^ + J{z^-h^)(l-h^) 



I{z) 



{l-z^Wz^-lP 



log- 



z + a/z^ — 6^ 



+ 



iog(i + vr^^) 1 , i + vr^&^ iog2 

log- 



and 



u:{z) 



1 



/X(l - ^2 

+ 



1 - Z 



log 



1 + ^ 



2(1-62) 



\-z 



z"" + h'' + ^{z^ -h^){l-lP) 
z^-b^ + ^{z^ -h^){l-h^) 



1-62 



log(i-62)-^ z 



'^2 _ 52 

1-62 



Since 

y{z) is simply given by 



/i = -log(l-62), 
W\z) - 2i2uj{z) 



1-^2 



log- 



Vz^ - 62 + 



Going back to the original u plane by the relation 



tanh 



u 



1 



e" + 1' 



the differential y{z)dz becomes 



l + e- + v/(l + e-")2-4e^-" 
y[z)dz = 2ctM log \ 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



Thus one recovers the equation given in jU] describing the mirror of a resolved conifold ^37j, 
which can be made identical to the equation of a deformed conifold by a change of variables. 
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5.3 The Special Geometry Relation and Free Energies 



The special geometry relation for this matrix model can be similarly obtained, although 
the branch cuts of W'{z), which run from —1 to — oo and from +1 to +oo, add an extra 
complication; one must specify the relative location of the branch cut of log|A| (taking 
Aoi = 0) in ()3.9|) and define the contour Bi accordingly. Then one differentiates Fq by /i and 
employ an identity similarly derived by using Riemann's bilinear relation. Although tedious, 
the calculation can be done similarly as before and we only describe the result: 

dFa r ri-e 



lim 



dzy{z) - W{1 - e) 



dyjz) 



dz 



(5.15) 
(5.16) 



where, as before, uj^j^{z) {uj^_^^_^^{z)) is a zero A-period abelian differential of the third 
kind having two simple poles, one at z = 1 {z = —1) with residue 1 and the other at z = 1 
{z = —1) with residue —1. Differentiating (j5.15j) by jj, again and using (j5.16|) . we find 



dij? 



Integrating over /i yields 
Fo = 



lim 



log 



1-e 



dz\ 



I -IP 



z + l 



log2e 



4(1-62)- 



E 

n=l 



H log 4 + const X /i + const. 



6 



(5.17) 



(5.18) 



ocs 



0-12 



agrees with the known result 



The genus-1 free energy may also be computed by the formula of Ambj0rn et. al. ^2] again: 



After some calculations we find 



- — (logMi+logJi + 41og26). 



Ml = Ji 



and hence obtain 



Fi = (^/i-log/i + log(l-e-'^)+41og2) 

roduces the ^ 
including the log term [HHl IBH] ■ 



(5.19) 



(5.20) 



(5.21) 



which correctly reproduces the known genus-1 Chern-Simons free energy -Fics = -^i + i| 
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6 Conclusions 



In this paper we have shown that a hermitian 1-matrix model with the potential and 
a unitary 1-matrix model with the potential Wu{U) are nothing but different descriptions 
of the same model if they are related with each other by the equation ()2.12j) . The extra 
logarithmic term on the hermitian side has been further interpreted as an effect of the 
coupling to 2N vectors. This duality is particularly useful to compute lower genus free 
energies in unitary as well as other matrix models with the unitary matrix Haar measure, 
and we have demonstrated the use of it in various examples. It allowed us not only to rederive 
the known facts by simply using the familiar hermitian matrix model technologies, but also 
to establish some new results on matrix models and gauge theory. We have shown that the 
unitary matrix formulation of the M = 2 pure SU (2) gauge theory correctly reproduces the 
genus-1 amplitude, and also given the rigorous forms of the special geometry relations in 
unitary as well as the Chern-Simons matrix models. 

This duality is not only of practical but also of conceptual importance; it has revealed that 
unitary and hermitian 1-matrix models have no essential difference and can be mapped to 
each other via the transformation (j2.12j) . An immediate consequence of this fact is that 
the partition function of a unitary 1-matrix model is also a r function of the KP hierarchy. 
It also provides an explanation of why the critical behaviors of both unitary and (even- 
potential) hermitian 1-matrix models are described by the same (Painleve II) equation when 
two clusters of eigenvalues join: this is accounted for by the universality with respect to even 
perturbations (e.g. j40ij) in hermitian 1-matrix models. 

It is amusing that while the topological B-model on a deformed conifold is known to be 
equivalent to c = 1 strings at the self-dual radius jH] and realized as the Penner model 02], 
the hermitian matrix model obtained after applying the duality to the Chern-Simons matrix 
model may be regarded as a "two-flavor" extension of the (generalized) Penner model. 

It would be interesting to apply this duality to the matrix models proposed to describe some 
five-dimensional gauge theories and compare the lower genus amplitudes with Nekrasov's 
formula 1^. 
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Appendix A The Multi-cut Free Energy 



In this Appendix we give a derivation of the formula (|3.9|) for the genus-0 multi-cut free 
energy in hermitian 1-matrix models. We start from the well-known expression |021 



,2 



— oo J — oo 



dXp{X)W{X)+ dX dX'p{X)p{X')\og\X~X'\ 



1 

fX"^ jJi J — OO 

which is valid for any number of the cuts n. It is convenient to define 



a{X) 



dX'p{X' 



Due to ()3.10|) we have 



a{ai) 
a{a2) 
' ' ' ) 

a{a2i) 



0, 



a{a2 



Cr{a2n) = 1- 

Using cr(A), the second term of ()A.1|) can be integrated by parts to give 

/oo " ^ ra2j 

dXp{X)Y, / (iA'p(A') log |A — A'l 

= / dXp{X)Y.[ a(A')log|A-A" ■ ' -'^^^^ 

J —oo • 1 \ . 



+ 



dX' 



a2j-i -"12^-1 

i-1 



A- A' 



(A.l) 



(A.2) 



(A.3) 



roo " / 1 1 

/ dXp{X) ^ - ^/ifclog |A - a2j\ Y /"fclog 1-^ - '^2j~i 

7 = 1 k=l k=l 



a2j 



-J2 dXa{X)-W'iX), 

7=1 -^aaj-i ^/i 



(A.4) 
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where in the last hne we have used the large N saddle point equation 

P(A') 



--w'i\)+2 r 



/i J —oo A — A' 

The last term is integrated by parts again to yield 



0. 



(A.5) 



E/ dXa{X)-W'{X) = - dXpiX)ViX)--(^{W{a,)-W{a,)) + 



Therefore we get 



Fo 



1 



2fi 



dXp{X)W{X) 



+ - (-^iW{a2) - W{a,)) + r p{X) log 

ji \ 2p J~oo A — 03 

+ ■■■ 

(-^(^K) - W^K+i)) + /_^P(A) loj 

+ ■■■ 

+ (^-^W^(a2n) + /^p(A)log|A-a2n|^ . 



A — a2j 



X- a2 



(A.7) 



Here the saddle point equation ()A.5|) can be used to replace a2k-i, a.2k with an arbitrary Aqa: 
on the fcth cut, for all /c = 1, . . . , n. With this replacement we end up with the formula (|3.9p . 



Appendix B A Derivation of Eg. (13. 181 ) 



Let S be a hyper-elliptic Riemann surface defined by 



2n 

Y[{z-ak) 

k=l 



(B.i: 



The first step in using Riemann's bilinear relation is to cut out S along the cycle 

An-iBn-iA-\B-\ ■ ■ ■ AiB.A^'B^' Br,P:'A^p-'B-'P£'A^P^' (B.2) 

to get a simply-connected polygonal region, which we denote by S'. For the definitions of 
the cycles see Figures 1 and 2. Next we consider 



tOi (i = 1, 



.n-1) 



(B.3) 



zo 
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Figure 2: More definitions of tlie contours. Tlie bold lines from ±i to Aq denote the branch 
cut of log(z T i)- ^605 are the same as those without tildes but lying on the second 
sheet. 



for some reference point zq on S', then Qi{z)^s are single- valued, everywhere holomorphic 
functions inside S'. Then for each i Riemann's bilinear relation can be obtained by computing 



9E' 



dzy{z)Qi{z) 



(B.4) 



in two ways, by the residue computation inside S', and by evaluating the integral along the 
boundary of S'. The total residue is trivially zero. On the other hand, one can evaluate the 
integral on the edges pairwise as 



+ , 

A, Ja-' 



+ , 



hm / + 



lim / + 



dzy{z)Qi{z 

dzy{z)Qi{z 
dzy{z)Qi{z 
dzy{z)Qi{z 

dzy{z)Qi{z 



0, 



n „ 

- V i W{z)uJi{z) - 47ri/ifii(-cx)), 

n „ 

f W{z)uji{z) + 47ri/ii7j(^oo), 
,=1 "'^j 



(B.5) 
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and 



dzy{z)Qi{z) 



dzy{z)Qi{z) = 0. 



Stokes' theorem implies that the sum of all these integrals is zero. Noting that 



lim 



we have 



— 47rz yU,- lim / tUi+f dzy(z) — f W(z)uJi(z) 



for i = 1, . . . , n — 1. 



(B.6) 



(B.7) 



(B.8) 



To derive a relation involving the period of a>„ we also consider 

/ dzy{z)Qn{z), 

JdT.' 

where 



(B.9) 



(B.IO) 



is also a single- valued holomorphic function inside S'. Taking special care about the singular 
behavior of y[z)VLn{z) near z = ±i, one may similarly compute (jB.9|l in two ways to obtain 



1 ^ r / " 

--J2f W{z)un{z) + lim -Tci ^ fij 



UJr, 



lim 

Ao—r — OO 



AJb. 



dzy{z) 



+ dziy{z) - W'iz))logiz' + 1) - ^WiAo) . (B.ll) 



One can then use ()B.8|) in this equation to see also that 



1 " 

--V^ W{z]LUi(Z 



lim 

Aq— > — CXD 



4 

1 



dzy{z) 



(B.12) 



for z = 1, . . . , n — 1, and hence, due to ()B.11|) . -i = 1, . . . , n. Since the last line is equal to 
1 



lim 



^ lim ( r f ' dzyiz) - W(t - e) - W(-i - e) 
4 Ao^-oo Vjao JAo 



B.13) 



we obtain the desired equation ()3.18p . 
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